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Abstract 

Detailed semi-empirical interaction potential calculations are performed to determine the potential 
energy surface experienced by the molecules CO, HCN, H2O and H2CO, when adsorbed on the 
basal plane (0001) of graphite at low temperature. The potential energy surface is used to find the 
equilibrium site and configuration of a molecule on the surface and its corresponding adsorption 
energy. The diffusion constant associated with molecular surface diffusion is calculated for each 
molecule. 

1 Introduction 

Observations in the cold regions of the interstellar medium (ISM) with typical temperatures of about 10 
- 20 K have revealed the presence of clouds which contain high densities (more than 10 3 particle per cm 3 ) 
of atoms, radicals and simple molecules such as H2, CO, CO2, HCN, H2O (the most abundant molecule), 
NH3, H2CO, CH4, CH3OH, as well as small solids typically composed of amorphous H2O, silicates and 
carbon grains in the form of graphite, amorphous structures and polycyclic aromatic hydrocarbon (PAH) 
molecules [H Hi- 
lt is nowadays thought that the dominant mechanism for the formation of molecules is provided by 
surface-catalysed chemical reactions on low temperatures interstellar dust grains. For instance, Watanabe 
and Kouchi [3] have experimentally produced the formaldehyde H2CO and methanol CH3OH molecules 
by successive hydrogenation of CO on the surface of icy grains at 10 K, as was theoretically proposed by 
Tielens et al. [4]. Recent laboratory experiments reported that more complex molecules such as glycine 
and other amino acids (which play an important role in the origin and evolution of life) as well as other 
organic species were formed from ultraviolet (UV) irradiation of the analogues of icy interstellar grains at 
low temperatures (T ~ 15 K) and containing small amounts of HCN, NH3, H2CO, and CH3OH species 

Theoretical quantum-chemical calculations on the synthesis of amino acids were reported by Woon 
[7] for UV-irradiated astrophysical ices and by Mendoza et al. [8] for ice analogues on the surface of a 
coronene (PAH) molecule. 

Investigations of the chemical processes leading to the production of amino acids on the surface of 
grains require a knowledge of the interactions between the molecules participating in these processes 
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on the surface, and also of their adsorption characteristics (such as the equilibrium configurations and 
diffusion constants) when they are isolated on the grain surface. 

The interaction of a single water molecule with a graphite sheet and with PAH molecules like the 
fused-benzene molecules (fbz)„ (C^Hq, C24H12, C54H18, C96H24,... for n = 1, 7, 19, 37,..., respectively) 
has been theoretically studied using ab-initio methods to determine the binding energies and equilibrium 
configurations [HI QH [HI HH EH [331 HSl HE] • For H 2 adsorbed on a graphite sheet, the calculated binding 
energies by various authors are -187 [TU], -72 [S], and -109 meV [17]. The available experimental value is 
-156 meV [IS]. 

For the water-benzene complex, many theoretical ab-initio calculations were reported. The calculated 
binding energies are around -120 [H] [TBI [14] , -169 meV [TTJ, and -99 meV [13]. The available experimental 
values are -98 [Uj, and -106 meV [20]. 

Feller and Jordan [12] have reported for water- (fbz) n complex up to n = 37 an energy binding -251 
meV which is twice as large as the recently reported value -126 meV by Lin et al. [13] . 

In the present work we study the adsorption properties of several of the molecules which are associated 
with the synthesis of amino acids. We use the graphite (0001) surface as a model substrate to represent 
the surface of an interstellar dust grain and we consider the surface adsorption of the various molecules 
on the model surface. We use an approach which has already been developed to calculate the infrared 
spectra of the NH3 ammonia molecule adsorbed on a variety of surfaces [211 [221 [23J [24] . The basic 
procedure involves the calculation of the potential energy surface for the adsorbed molecule as it moves 
on the adsorbing surface. The environmental effect on the inversion spectrum of NH3, for example, can 
be simulated by changes in the depth or curvature of an effective double well potential. Depending on the 
substrate surface it may happen that, for the NH 3 molecule, the central barrier height increases (giving 
reduced splittings within each of the lowest two vibrational levels) or an asymmetry appears between 
the two potential wells. Even a slight asymmetry can inhibit the inversion motion in NH3, so that the 
eigenstates are effectively localized on one side or the other of the central potential barrier. These earlier 
calculations with NH3 have given us experience of the calculational techniques involved and they led to 
calculated spectra which were in agreement with the experimental ones. 

The isosteric heat of adsorption for NH3 adsorbed on the graphite substrate calculated in our model 
was 178 meV [24], in agreement with the experimental value of 172 meV given by Avgul and Kiselev [18] 
at a temperature T — 195 K. 

The various parts of section 2 of this work describe the important contributions to the potential energy 
of the adsorbed molecule on the adsorbing surface. Section 3 describes the minimization process needed 
to find the possible equilibrium positions of the adsorbed molecule on the surface and gives details of 
the observables which are to be calculated. Section 4 gives the results obtained for the several molecules 
studied; typical potential surfaces are shown in the figures, to augment the sometimes complicated verbal 
description of the equilibrium positions of the molecule. 

Two tables set out some important calculated numerical parameters for the various adsorbed molecules. 
An appendix sets out some of the mathematical transformations needed to transform various quantities 
between the internal molecular axes and the chosen reference set of axes needed to give a unified descrip- 
tion of the various quantities. 

2 Theoretical background 

An adsorbed species can be either physically or chemically bonded to the surface of the substrate crystal. 
In physisorption the interaction is due to the comparatively weak van der Waals-London and electric field- 
multipole forces, and the interaction potential energy is fairly low (< 0.5 eV), with no charge exchange 
between the adsorbed species and the substrate. 

2.1 The physisorbed molecule on a graphite substrate 

In this work we consider a single molecule such as CO, HCN, H2O or H2CO physisorbed on a graphite 
basal plane (0001) at low temperatures. 

The graphite crystal is extensively used as a model to simulate the behaviour of interstellar grains in 
studies of the physical and chemical behaviour of adsorbed molecular species. Experimentally detected 



2 



absorption features due to carbon grains have been attributed to the presence of small spherical graphite 
crystals of size ~ 200 A in interstellar dust clouds [TJ [5] • 

In order to reduce the required numerical calculation times, the graphite crystal is modelled by using 
ten planes with radii ~ 35 A (~ 1500 carbon atom per plane). This does not significantly influence 
the interaction potential energy between the molecule and the substrate (giving an error of less than 2 
percent). 

Figure 1 shows the geometrical characteristics of a molecule (H 2 CO for instance) adsorbed on the 
graphite surface. 

To perform the calculations, we need to define: 

i) an absolute frame (0,X,Y,Z) connected to the surface of the substrate, with respect to which 
all of the vectors and the tensors associated with the molecule and the graphite carbon atoms must be 
described; 

ii) a frame (G,x,y,z) tied to the molecule with centre of mass G, and with respect to which the 
molecular vibrational motions are still described. The rotational matrix transformation from the absolute 
frame to the molecular frame is given in Appendix A, in terms of the molecular orientational degrees of 
freedom Q, = (ip,9,x) (precession, nutation, and proper rotation Euler angles). 



2.2 Interaction potential energy 

The interaction potential energy Vmg between an adsorbed molecule and the graphite substrate can be 
written as the sum of three parts 

V MG = V LJ + V E + V l , (1) 

which are detailed below. 



2.2.1 Quantum contributions 

The first term Vlj in Eq. {T]) is generally treated by using semi-empirical 12-6 Lennard-Jones (LJ) pairwisc 
atom-atom potentials to describe the short-range atomic repulsion energy arising from the overlap of the 
electronic clouds of the adsorbate and the substrate at very short distances and also the long-range many- 
body dispersion energy arising from the lowering of the quantum zero-point energy of the system by the 
correlated fluctuations of the adsorbate and substrate atomic dipoles. 

However, it is widely known that carbon atoms in the graphite basal plane (0001) exhibit anisotropic 
polarizability effects which modify the spherical shape of these atoms. To account for these anisotropies, 
Carlos and Cole [25] proposed a modified expression for Vlj , of the form 
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1 + 7a(1 - -cos 2 %) 



(2) 



In this equation, r.y is the distance vector between the ith atom of the molecule and the jth carbon 
atom of the graphite substrate (Figure 1), while and aij are the mixed LJ potential parameters, 
obtained from the usual Lorentz-Berthelot combination rules = sjeaejj and 2ct 2 j = era + ajj . The 
two coefficients 7^ and ja specify the modifications to the repulsion and dispersion parts of the energy, 
respectively, and 9ij is the angle between the surface normal axis and the r^- vector. The various molecule- 



graphite mixed parameters are given in Table 1. 
been improved several times in our other works. 



Note that the LJ parameters used in this work have 



2.2.2 Electrostatic and induction contributions 

It is well known that the bonds in the graphite crystal give rise to aspherical atomic charge distributions 
which produce multipole electrostatic fields with rapid spatial variations external to the graphite [9]. 
This phenomenon is modelled by introducing axially symmetric quadrupole moment tensors O c at each 
carbon atom site. 
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Figure 1: Geometrical characteristics of an H2CO molecule adsorbed on the graphite (0001) substrate. 
(0,X,Y,Z) and (G,x,y,z) represent the surface absolute frame and the molecular frame, respectively. L 
and a are the interlayer spacing and the distance between two next-nearest neighbour carbon atoms of 
the graphite and qf , q| and f3 e define the equilibrium internal coordinates of the adsorbed molecule. 
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Table 1: Lennard- Jones parameters for atom- graphite combinations. 



Gr - Atom 


e (meV) 


a (A) 


(a) 

7a 


(a) 

7r 


Gr - H 


2.265 


2.965 


0.4 


-0.54 


Gr - 


3.450 


3.141 


0.4 


-1.05 


Gr - C 


2.981 


3.305 


0.4 


-1.05 


Gr - N 


2.811 


3.390 


0.4 


-1.05 



(a) From Refs. [2SJ [H [23 HO EU . 



The quadrupole moment tensor Q c of the jth carbon atom produces the electrical potential 

^» = iL p xw(^):0 C) (3) 

on the molecule at its centre of mass position r. In turn, the multipole moments (/x, S,...) of the molecule 
produce the electrical potential 

*"(r,) = L p x V (^) .„ + \ Lp x VV (^) : + (4) 

on the jth carbon atom at position Tj . 

In Eqs. ([3]) and Q, L p is a screening factor, equal to 1 for the surface plane and to 2/(e + 1) for the 
internal planes (p > 2), where e is the static dielectric constant of the graphite crystal, which we take to 
be 2.8. 

The electrostatic and induction contributions to the interaction potential energy Vmg are 

V E = V tz.V^r) + 1 V 6 : VV$ J (r) + ... (5) 

O 

J 3 

* = -\Y. [ V * M fe) : « C : V$M fe)] - \T. [V^(r) : « : V^'(r)] , (6) 

where a and a c are the polarizability tensors of the molecule and of the graphite carbon atoms, respec- 
tively. 

The elements of the polarizability a c and quadrupole moment C tensors (Eqs. (J3HHJ) ) of the carbon 
atoms are given in the absolute frame (0,X,Y,Z) by 

a xx = a YY = a l and "zz = «|j i ( 7 ) 

©xx = e YY = -ie c zz = -ie c , (8) 

with all other elements being zero. 

Moreover, it must be noticed that the molecular frame corresponds, for each of the studied molecules, 
to its principal frame of inertia, in which only the z-component of the dipole moment /i vector and the 
three diagonal elements of the polarizability a and quadrupole moment O tensors are non-vanishing 
elements. In the potential energy calculations these quantities must be expressed in the absolute frame 
by using a rotational matrix transformation (Appendix A). 

The various geometrical and electrical parameters of the graphite substrate and of the studied 
molecules are given in Table 2. 

2.2.3 Separation of Vmg 

The distance vector in Eq. ([2]) can be expressed as 

r y = T i _r _ Vu ( 9 ) 
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Table 2: Molecular and graphite parameters: internal bonds and angles, dipole moments, polarizabilities, 
quadrupole moments, and rotational constants for the molecules studied. 



Molecule 


CO 


HCN 


H 2 


H 2 CO 


Graphite^ 




{q e } (A) 


1.144 


1.150 


0.958 


1.198 


L (A) 


3.36 




- 


1.079 


_ 


1.121 


a (A) 


2.46 


P c (deg) 


_ 


_ 


54.7 


58.1 






fi (D) 


0.112 


2.986 


1.855 


2.330 






«xx (A 3 ) 


1.63 


2.05 


1.53 


2.51 


al (A 3 ) 


1.44 


"yy (A 3 ) 


1.63 


2.05 


1.42 


1.94 


a|| (A 3 ) 


0.41 


a zz (A 3 ) 


2.60 


3.66 


1.47 


2.90 






e xx (DA) 


1.02 


-2.20 


2.63 


1.91 


e c (dA) 


1.0 


©yy (DA) 


1.02 


-2.20 


-2.50 


-5.01 






e zz (dA) 

A (cm -1 ) 


-2.05 


4.40 


-0.13 


3.09 






1.93 


1.48 


27.33 


1.22 






B (cm- 1 ) 


1.93 


1.48 


14.57 


1.13 






C (cm" 1 ) 






9.49 


9.40 







(a) From Refs. [25j HE W\ H3 EU • 

where r and Tj are the position vectors of the molecular centre of mass and the jth graphite 
carbon atom, respectively. 

The vectors rji characterize the positions of the atoms in the molecule. They, as well as 
the molecular electrical parameters fi, 9, a, depend on the molecular internal vibrational 
motions. They could be written, with respect to the (G,x,y,z) frame, as 

m = vt + ^T J a»Q v + 

V 

e = o e + J2^Q» + -, 

v 

a = a e + J2<l"Q» + -- (10) 

In these expressions the superscript e refers to the equilibrium internal configuration of the 
molecule (rigid molecule) and a", h" , c v , and d" are the first derivatives of rji, fj,, 9, and 
a with respect to the normal coordinate Q u associated with the vt\i molecular vibrational 
mode with frequency lo u . 

However, the present work is devoted to finding the adsorption energies and diffusion 
constants of single molecules adsorbed on the graphite substrate. Therefore, we can assume: 
i) the substrate to be rigid ({r^} = {r|}), H) the adiabatic approximation to be valid, 
permitting us to separate the high frequency vibrational modes of the molecule from its low 
frequency external translational and orientational modes, and Hi) the dynamical coupling 
between all of the molecular and the graphite degrees of freedom to be negligible. 

The interaction potential energy Vm.g can then be written as 

V MG - V MG (r, fl, {Ql}) + V mG ({Q v }, r e , fl e ), (11) 

where Vmg (r, {Qt}) represents the external motion-dependent part of the potential energy 
experienced by the non-vibrating molecule. Vmg({<3i/}, r e , Sl e ) characterizes the vibrational 
dependent part, for the molecule at its equilibrium position and orientation; in a perturba- 
tive approach it contributes to the molecular gas phase vibrational Hamiltonian and leads 
to small vibrational frequency shifts. 
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3 Adsorption observables 



At sufficiently low temperatures and within the rigid molecule and substrate approximation, 
the potential energy surface Fmg(-^,^)i which is experienced by the adsorbed molecule as 
it moves laterally above the surface can be calculated by minimizing Vmg with respect to 
both the perpendicular distance Z between the molecular centre of mass and the surface 
and the molecular angular coordinates (ip,0,x)- 

The resulting energy map gives information about the equilibrium site (r e ,f2 e ), and the 
lateral diffusion valley leading from one equilibrium site to an adjacent one. This makes 
it possible to determine the isosteric heat of adsorption or adsorption energy and also the 
surface diffusion constant. 



3.1 Adsorption energy 

The isosteric heat of adsorption E a of a single adsorbed molecule on the graphite substrate 
can be defined as the energy required to keep the molecule with thermal energy (~ + l)kT at 
its equilibrium configuration on the surface of the substrate. It is approximately expressed 
as 




where n is the number of external degrees of freedom for the molecule (5 for linear molecules 
and 6 for non linear molecules) and k is the Boltzmann constant. The w s (s=l to n) are 
the frequencies associated with the translational and orientational motions of the molecule 
around their equilibrium positions. 



3.2 Surface diffusion constant 

Let Vmg(t, £,~){t)) be the instantaneous potential energy surface experienced by the adsorbed 
molecule with mass m as it moves along the diffusion valley described by the coordinate path 
r. £ 7 (t) describe the instantaneous distance from the surface (7 = Z) and the instantaneous 
orientation (7 = tp, 6, x) of the molecule for the given r value. 

If we disregard the dynamics of the £ degrees of freedom around their r-dependent equi- 
librium values £ along the diffusion motion of the molecule and use the classical transition 
state theory |30j, then the jumping rate for the diffusion from one equilibrium site Si to an 
equivalent adjacent one S2 can be written as |31| 

K (T) ~ ( kT ) * exphAy MG (r )/fcr] 



\ 2nm* 



J dTexp[-AV MG (T)/kT] 



where AVmg(t) = Vmg(t) Vmg(ti), and To corresponds to the saddle point position along 
the diffusion valley. The diffusion constant is then 

D{T) = X 2 K(T), (14) 

where A = T2 — T\ is the jump length, i.e., the distance (in the diffusion direction) between 
the two adjacent equilibrium sites. 

In Eq. (|13p m* is an effective molecular mass which is due to an inertial effect resulting 
from the changes in the distance of the molecule from the surface and in the molecules 
orientation as the diffusion process proceeds. It is expressed as 

m* = m(l + Y\—al\, (15) 
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where A 7 represents the molecular mass m for 7 = Z and the molecular moments of inertia I a , 

lb, and I c for 7 = ip, 0, and x, respectively. The a 7 = ^f- represent deformation parameters 
as the molecule migrates from an equilibrium site to the saddle point position along the 
diffusion valley. 

4 Numerical results 

4.1 Potential energy surfaces and equilibrium configurations 

The potential energy surfaces Vmg(X,Y) for the CO, HCN, H2O, H2CO single molecules 
adsorbed on the graphite substrate are presented in Figures 2 using a square surface area 
of (a x -^|), where the lattice parameter a is shown in Figure 1. 

In Table 3 we give, for each adsorbed species, some characteristics of the adsorption, 
such as the most favourable positions and orientations and their associated energy minima 
^mc as wen as tne energy maxima V^ G . 

4.1.1 CO admolecule 

The potential energy surface experienced by the CO adsorbed molecule is presented in 
Figure 2a. The molecule exhibits an equilibrium configuration with an associated energy 
minimum V^ G = -120.2 meV located above the hexagon site centre at a distance Z — 3.08 
A from the surface and for a nearly flat orientation (0 e ~ 95 deg) position. 

The perpendicular orientational motion around this equilibrium configuration is a li- 
brational motion with an energy barrier height of about 40 meV, while that in the plane 
parallel to the surface (cp) is a free rotational motion. 

The energy maximum V^ G = -106.5 meV is obtained with the molecule in a nearly fiat 
configuration at a distance Z — 3.24 A above the graphite carbon sites. There is a surface 
corrugation energy V™ G — V^ G = 13.7 meV. 

The saddle point of the diffusion path is reached when the molecular centre of mass is 
above the middle of the carbon-carbon bonds at Z — 3.20 A from the surface, with the 
molecule perpendicular to the C-C bonds and being always in the nearly flat configuration 
(Figure 3a). The associated energy barrier height is AVmg(to) = 12 meV (Table 6). 

4.1.2 HCN admolecule 

The linear HCN adsorbed molecule presents a potential energy surface with six most 
favourable sites per graphite hexagon (Figure 2b). The energy minimum V^ G — -175.5 
meV is reached with the molecule in a nearly flat configuration (9 e ~ 100 deg) with its 
centre of mass at Z — 3.20 A on the half of the C-C bonds, with a displacement of 0.28 
A inside the hexagon surface and with the hydrogen atom of the molecule pointing in the 
hexagon site centre direction. 

In these equilibrium sites there is a strongly hindered orientational motion perpendicular 
to the surface (9 motion) with an energy barrier height of about 75 meV and a moderately 
hindered parallel one ((p motion) with an energy barrier height of about 16 meV. 

However, it is interesting to note that this latter orientational motion, combined with 
the translation motion of the molecule between the six equilibrium sites along a quasi- 
circular trajectory of radius ~ 0.95 A around the hexagon site centre, could be regarded 
as a quasi-free motion (Figure 3b) with an energy barrier height < 1.5 meV. Thus, an 
associated effective moment of inertia must be defined with respect to a new HCN y-axis 
lying at 0.95 A from its centre of mass (between the C and H atoms). 

The energy maximum of V^ G — -162.5 meV is obtained with the molecular centre of 
mass at Z — 3.24 A, nearly above the hexagon site centre and also with the nearly flat 
configuration. This induces a surface corrugation energy of 13 meV. 
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Table 3: Adsorption characteristics for CO, HCN, H2O and H2CO molecules. 



Molecule 


X° (A) 


y° (A) 


Z c (A) 


tp c (deg) 


C (deg) 


X° (deg) 




V^ G (mcV) 


CO 


0. 


0. 


3.08 


free 


95 




- 120.0 


- 106.5 


HCN 


- 0.95 


0. 


3.20 





103 




- 175.5 


- 162.5 


H 2 


- 1.23 


- 0.36 


3.00 


90 


100 


90 


- 138.0 


- 129.5 


H 2 CO 


- 1.11 


0. 


3.04 





90 


90 


- 208.5 


- 186.5 



Moreover, this position also corresponds to the saddle point of the diffusion motion from 
one equilibrium site of a hexagon to one equivalent equilibrium site of an adjacent hexagon 
(Figure 3b). 

4.1.3 H2O admolecule 

The potential energy surface of the non-linear H2O molecule adsorbed on the graphite 
substrate is presented in figure 2c. The most favourable sites are obtained for the molecular 
centre of mass at Z — 3.00 A above the C-C bonds, with a displacement of 0.36 A from 
the C graphite carbon atom, the oxygen atom of the molecule being always close to this C 
graphite carbon atom (Figure 3c). 

The energy minimum V^ G — -138.0 meV is obtained with the molecule in a nearly fiat 
configuration (9 e ~ 100 deg and \ — §) an d sitting astride the C-C bonds (Table 3), and 
is in agreement with the available experimental value -156 meV obtained by Avgul and 
Kiselev [18] . 

It should be noted that in the equilibrium sites the angular motions of the C2 molec- 
ular symmetry axis (z-axis) parallel (ip motion) and perpendicular (9 motion) to the sur- 
face around their equilibrium values are moderately hindered motions with energy barrier 
heights of about 8 and 15 meV, respectively. The angular (spinning) motion about this 
axis (\ motion) is a completely hindered motion, the energy barrier height being about 400 
meV. 

The energy maximum V^ G — -129.5 meV is reached for Z — 2.92 A above the hexagon 
site centre with the C2 molecular symmetry axis perpendicular to the surface of the graphite 
substrate. The surface corrugation energy is then 8.5 meV. 

In Figure 3c we sketch the valley of diffusion of the water molecule on the surface. 
The barrier height at the saddle point is AVmg(to) — 1 meV (Table 6). Note however, 
that to obtain such a motion it is necessary to combine the translational motion and a 
reorientational motion parallel to the surface (ip — ±-f )• 

4.1.4 H 2 CO admolecule 

Like the HCN molecule, the H2CO molecule adsorbed on the graphite substrate exhibits 
six equilibrium sites per hexagon (Figure 2d). The energy minimum V^ G — -208.5 meV 
is obtained for the molecule in a flat configuration (8 e — 90 deg and \ e — 90 deg) with 
its centre of mass at Z — 3.04 A on half of the C-C bonds, with a displacement of 0.12 A 
inside of the hexagon surface and with its C2 symmetry axis (z-axis) pointing towards the 
hexagon site centre. 

It should be noted that in its equilibrium sites H2CO behaves nearly like HCN with 
regard to the perpendicular 9 and parallel ip angular motions of its z-axis, in one hand, but 
like H2O with regard to its spinning motion about this axis. 

The energy maximum is V^g — -186.5 meV, with the molecule in the flat configuration 
at Z — 3.12 A above the hexagon site centre. The surface corrugation energy is then 22 
meV. 

The saddle point of the diffusion path is obtained for the molecule in the flat configuration 
at Z — 3.04 A above the C carbon site with a corresponding energy barrier height AVmg(to) 
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Table 4: Calculated frequencies (cm -1 ) associated with the translational and orientational oscillations of 
CO, HCN, H 2 and H 2 CO molecules. 



Molecule 




LUy 






UJ0 




CO 


26.6 


26.6 


62.1 


~ 


76.6 




HCN 


21.0 


23.4 


79.0 


32.3 


109.7 




H 2 


~ 


~ 


90.3 


22.6 


35.5 


233.9 


H 2 CO 


25.0 


~ 


155.7 


21.0 


91.1 


203.3 



Table 5: Calculated adsorption energy E a (meV) of CO, HCN, H 2 and H 2 CO molecules for three typical 
temperatures. 



Molecule 


T = 10 K 


T = 40 K 


T = 100 K 


CO 


111.3 


114.8 


114.4 


HCN 


168.5 


163.8 


159.5 


H 2 


117.6 


123.7 


129.6 


H 2 CO 


181.0 


186.6 


189.8 



= 5 meV (Table 6). Moreover, the rotation of the molecule parallel to the surface plane is 
a free motion on this site. 

We note that the fact that the equilibrium configuration is nearly above the carbon- 
carbon bonds for HCN, H 2 0, H 2 CO molecules results from their strong multipole moments. 
Their quantum V^j, induction Vi and electrostatic Ve contribution parts represent ~ 80%, 
~ 15% and ~ 5%, respectively, of the total potential energy for these molecules. 

The weakly polar CO molecule behaves as a quadrupolar molecule and adopts an equi- 
librium configuration above the hexagon site centre. Its quantum contribution part Vlj 
represents more than ~ 95% of the total potential energy. 

4.2 Adsorption energy calculations 

The calculation, using Eq. (1121) . of the adsorption energy values for single molecules ad- 
sorbed on the graphite substrate requires a knowledge of the frequencies oj a associated 
with the linear and angular oscillations associated with the translational and orientational 
motions of the molecules around their equilibrium configurations on the most favourable 
sites. 

A discrete variable representation method |32j was used to solve the Schrodinger equa- 
tion giving the energy levels of a molecule undergoing linear and angular oscillations at 
its equilibrium site, with the other degrees of freedom being assumed to remain at their 
equilibrium values. The various frequencies calculated for the molecules treated in this 
work are reported in table 4. 

As can be expected from the potential energy surfaces, the frequencies associated with 
the orientational motions perpendicular to the surface (8 and \ motions) (except the 9 
motion for H 2 molecule) are large compared to those associated with the orientational 
motion parallel to the surface (ip motion). 

Finally, the calculated frequencies were introduced into Eq. (|12[) and the adsorption 
energies were determined for three typical temperatures. The resulting values are reported 
in table 5. 
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Table 6: Calculated surface diffusion constant D (cm 2 /s) of CO, HCN, H2O and H2CO molecules for 
three typical temperatures. 



Molecule 


in 
m 


AVmgM (meV) 


T = 10 K 


T = 40 K 


T = 100 K 


CO 


1.01 


12. 


3.3 xl0~ lu 


1. xlO" 5 


7.3 xl0~ b 


HCN 


1.002 


13. 


1. xl0~ 10 


1.2 xl0~ 5 


7.2 xl0~ 5 


H 2 


1.11 


1. 


1.5 xl0~ 5 


4.3 xl0~ 5 


7.2 xl0~ 5 


H 2 CO 


1.27 


5. 


2.2 xl0~ 7 


1.8 xlO" 5 


5.1 xl0~ 5 



4.3 Diffusion constant calculations 

The study of the potential energy surface for each adsorbed species allows us to determine 
the diffusion path, its length A and potential energy function AVmg(t), and the deformation 
parameters a 7 as the molecular diffusion proceeds on the surface. 

It must be mentioned that large values for a 7 correspond to important equilibrium 
changes when the molecular centre of mass migrates from the minimum to the saddle point 
position along the diffusion valley. For instance, the angular deformation parameter values 
a v are 1.03 and 0.74 rad/A for H 2 and H 2 CO admolecules, respectively. 

In contrast, the angular deformation parameters ag and a x are negligibly small for all 
the admolecules because of the very weak changes of their perpendicular orientation with 
respect to the surface during the diffusion motion. 

The linear deformation parameter values a z are 0.10, 0.05, 0.04 and 0.00 for CO, HCN, 
H 2 and H 2 CO admolecules, respectively. 

Finally, the possible diffusion paths are indicated in figures 3 and the corresponding 
energy barrier heights AVmg(to) are given in table 6. 

Try 

Introducing the calculated deformation parameters in Eq. (|15|) the mass ratios are 

TO 

obtained and presented in table 6. Note that for the CO and HCN admolecules the mass 
changes are negligibly small, while for H 2 and H 2 CO the effective mass increases by 11% 
and 27%, respectively. 

5 Discussion 

In the present work we have studied the adsorption of several molecules on a graphite 
substrate; the molecules were selected because of their biophysical interest, as explained in 
the introduction. 

The major calculation was that of the potential energy surface for each molecule and 
the figures show the resulting surfaces. A minimization calculation then gave the possible 
equilibrium positions for each molecule, while the potential energy surface indicated the 
favourable pathways between these positions which would facilitate surface diffusion of the 
adsorbed molecule. 

The calculated adsorption energies as shown in table 5 show only a relatively weak 
variation with temperature, although the particular values of the molecule-surface interac- 
tion parameters for the HCN molecule lead to a temperature dependence which is in the 
opposite direction to that for the other three molecules studied. 

The calculated diffusion constants for the molecules are shown in table 6 and it is clear 
that the different saddle point barrier heights and effective masses lead to a variation of 
several orders of magnitude between the low temperature diffusion constants of the various 
molecules. In effect, one can remark that at T — 10 K the H 2 molecule could diffuse 3 
orders of magnitude more easily than H 2 CO and 5 orders of magnitude more easily than 
CO and HCN molecules. 

However, as the temperature increases the thermal energy kT becomes sufficiently large 
to overcome the barriers for all the molecules, with the consequence that the diffusion 
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Figure 3: The calculated equilibrium positions for the adsorbed molecule on the graphite substrate, 
indicating the energetically favoured paths for diffusion between the equilibrium positions, a (CO), b 
(HCN), c (H 2 0) and d (H 2 CO). 
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constants vary much less from molecule to molecule as the temperature rises. 
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Appendix : Rotational matrix transformation 

The unitary matrix M characterizing the transformation from the surface absolute frame 
(0,X,Y,Z) into the molecular one (G,x,y,z) through the Euler angles tp, 9 and \ i s given by 

m 

cos</3 cos6> cosx - simp sin% sintp cos# cosx + cos<^ sin% - sin 9 cosx \ 

cos <p cos 9 sin \ ~ sm <P cos X ~ sm f cos $ sm X + cos f cos X sm $ sm X 

cos <p sin 9 sin ip sin 9 cos / 

References 

[I] Draine, B. T. Astrophys. J. 1988, 333, 848. 

[2] Mathis, J. S.; Whiffen, G. Astrophys. J. 1989, 341, 808. 

[3] Watanabe, INT.; Kouchi, A. Astrophys. J. Lett. 2002, 571, 173. 

[4] Tielens, A. G. G. M.; Whittet, D. C. B. in Molecules in Astrophysics; Ed.; van Dishoeck 
E. F. (Dordrecht: Kluwer) 1997, 45. 

[5] Bernstein, M. P.; Dworkin, J. P.; Sandford, S. A.; Cooper, G. W.; Allamandola, L. J. 
Nature 2002, 416, 401. 

[6] Munoz Caro, G. M.; Meierhenrich, U. J.; Schutte, W. A.; Barbier, B.; Arcones Segovia, 
A.; Rosenbauer, H.; Thiemann, W. H. P.; Brack, A.; Greenberg, J. M. Nature 2002, 416, 
403. 

[7] Woon, D. E. Astrophys. J. Lett. 2002, 571, 177. 

[8] Mendoza, C.; Ruette, F.; Martorell, G.; Rodriguez, L. S. Astrophys. J. Lett. 2004, 601, 
59. 

[9] Vernov, A.; Steele, W. A. Langmuir 1992, 8, 155. 

[10] Gale, R. L.; Beebe, R. A. J. Phys. Chem. 1964, 68, 555. 

[II] Feller, D. J. Phys. Chem. A 1999, 103, 7558. 

[12] Feller, D.; Jordan, K. D. J. Phys. Chem. A 2000, 104, 9971. 

[13] Lin, C. S.; Zhang, R. Q.; Lee, S. T.; Elstner, M.; Frauenheim, Th.; Wan, L. J. J. Phys. 
Chem. B 2005, 109, 14183. 

[14] Li, S.; Cooper, V. R.; Thohauser, T.; Puzder, A.; Langreth, D. C. J. Phys. Chem. A 
2008, 112, 9031. 

[15] Tsuzuki, S.; Honda, K.; Uchimaru, T.; Mikami, M.; Tanabe, K. J. Am. Chem. Soc. 
2000, 122, 11450. 

[16] Zimmerli, U.; Parrinello, M.; Koumoutsakos, P. J. Chem. Phys. 2004, 120, 2693. 



M(cp,e, x ) 



14 



[17] Karapetian, K.; Jordan, K. D. Water in Confined Environments; Devlin, J. P.; Buch, 
V., Eds.; Springer: New York, 2003, 139. 

[18] Avgul, N. N.; Kiselev, A. V. Chem. Phys. Carbon; Walker, P. L., Ed.; Dekker: New 
York, 1970, 6, 1. 

[19] Cheng, B.-M.; Grover, J. R.; Walters, E.A. Chem. Phys. Lett. 1995, 232, 364. 

[20] Courty, A.; Mons, M.; Dimicoli, I.; Piuzzi, F.; Gaigeot, M.-P.; Brenner, V.; Pujo, P. 
d.; Millie, P. J. Phys. Chem. A 1998, 102, 6590. 

[21] Lakhlifi, A.; Picaud, S.; Girardet, C.; Allouche, A. Chem. Phys. 1995, 201, 73. 

[22] Lakhlifi, A.; Girardet, C. J. Chem. Phys. 1996, 105, 2471. 

[23] Lakhlifi, A. Eur. Phys. J. D 2000, 8, 211. 

[24] Lakhlifi, A.; Killingbeck, J. P. J. Phys. Chem. B 2005, 109, 11322. 
[25] Carlos, W. E.; Cole, M. W. Surf. Sci. 1982, 119, 21. 
[26] Vidali, G.; Cole, M. W. Phys. Rev. B 1984, 29, 6736. 

[27] Hansen, F. Y.; Bruch, L. W.; Roosevelt, S. E. Phys. Rev. B 1992, 45, 11238. 

[28] Hansen, F. Y.; Bruch, L. W. Phys. Rev. B 1995, 51, 2515. 

[29] Ionov, S. I.; La Villa, M. E. J. Chem. Phys. 1992, 97, 9379. 

[30] Voter, A. F.; Doll, J. D. J. Chem. Phys. 1984, 80, 5832. 

[31] Lakhlifi, A.; Girardet, C. J. Chem. Phys. 1991, 94, 688. 

[32] Light, J. C; Hamilton, I. P.; Lill, J. V. J. Chem. Phys. 1985, 82, 1400. 

[33] Rose, M. E. Elementary Theory of Angular Momentum; Wiley: New York, 1967. 



15 



